NONSTANDARD FEYNMAN PATH INTEGRAL 
FOR THE HARMONIC OSCILLATOR 
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Abstract. Using Nonstandard Analysis, we will provide a rigorous computation 
for the harmonic oscillator Feynman path integral. The computation will be done 
without having prior knowledge of the classical path. We will see that properties of 
classical physics falls out naturally from a purely quantum mechanical point of view. 
We will assume that the reader is familiar with Nonstandard Analysis. 



I. Introduction. 

In quantum mechanics, we are interested in finding the wave function which 
satisfies Schrodinger's equation. Equivalently, we can find the propagator or integral 
kernel K(q, q , t) which satisfies 



. k dK(q,q ,t) 
in dt 



2m 



A 9 + V(q) 



K(q,qo,t), 



K(q,q ,0) = S(q-q o ), q,q eR d . 
Formally, the wave function is related to the propagator via 

-+oo 



ip(q,t) 



K(q,q ,t)(j)(q )dq = ( K(q,q ,t), (j)(q ) 



with boundary condition (p(q, 0) = (p(q). In Feynman's formulation of quantum 
mechanics, he proposed that the propagator is given by a functional integral, also 
referred to as Feynman path integral or just path integral in physics literature, 



;i.3) 



K(q,qo,t) 



:{t)=q 



exp 



a;(0)=qo 



lS[X h {s)] }dx(s) 



lim 

n — >oo 



w d ,n exp 



—S {x n+ i . . . xo} 
a 



dxi . . . dx n , where 



x = q ,x n+ i =g,e=-, 

n 



w d 



d(n+l) 

m \ — 2 — 



2inhe. 

n+l 

S {x n+1 ...x }=} j 

3=1 



m I Xj — Xj-i 
~2 



~V(xj) 
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The integrals in the second line of (1.3) are (/-dimensional improper Riemann in- 
tegrals. The integral in the first equality is the purely formal path integral which 
integrates over all paths x (s) , < s < t, with x (0) = qo, an d x(t) = q. The 
quantity S [x (s)] is the action integral 



(1.4) 



[x («)] = f 
Jo 



L (x, x) ds, 



m 



L (x, x) = —x — V(x) . 



The motivation for the notation of the formal path integral is that as n goes to infin- 
ity, the sum in the exponent becomes the action integral. In this manner, Feynman 
was able to deduce Classical Mechanics from Quantum Mechanics through the ac- 
tion integral and H — > 0. In this paper, we will compute the harmonic oscillator path 
integral without prior knowledge of classical physics. It turns out that properties 
of classical physics naturally falls out of the computation without taking h — > 0. 
The computation shows that the path integral separates into the product of two 
quantities, one independent of H, the other dependent on h. The quantity which is 
independent of h contains properties of classical physics. Thus, in some sense we 
are deviating from the standard interpretation of the path integral being a bridge 
between Quantum Mechanics and Classical Mechanics via the action integral and 
h. We are considering Quantum Mechanics as purely Quantum Mechanics and 
extracting properties of Classical Mechanics without prior knowledge of Classical 
Physics. 

Mathematically, the formal integration over paths can not be a rigorously well 
defined measure theoretic integration because of the oscillatory nature of the inte- 
grand(see [1] and [2]). A popular technique to make sense of the (1.3) is to replace 
t by —it and use the Wiener integral(see [3] and [4]). 

In Nonstandard Analysis, we have that lim,^^ a n = a iff * for any 

infinite natural number uj G *N — N, with {*a m } m6 * N being the *-extension of 



and ~ means that a w + h u 



a where h u an infinitesimal. 



We can use 



Nonstandard Analysis to define the path integral; the standard part can replace the 
limit in (1.3). Using Nonstandard Analysis to replace the limit in (1.3) is not a new 
concept (see [5], [6] and references within), doing so partially solves the problems of 
the Feynman path integral on the propagator. 

We can redefine (1.3) in the following manner: let u E *N,m,t G M+,e = 
*V(x) : *M. d — > *R be an internal function, and xq = qo 7 x u+ i = q be fixed points 
R d . We call the expression: 



t_ 
in 



(1.5) 



(- 



m 



\ 2niht 



J exp 



16 
Jl 



E 



m 
~2 



dxi . . . dx u 



an internal functional integral. In (1.5), all integrals are *-transformed improper 
Riemann integrals. In particular, if for all uj G *N — N, the standard part of the 
internal functional integral exists and it is independent of the choice of uj, we call 
the standard part a standard functional integral or Feynman path integral and 
denote it by 

(1.6) 



NONSTANDARD FEYNMAN PATH INTEGRAL FOR THE HARMONIC OSCILLATOR 3 



Equation (1.6) is just a Nonstandard Analysis way of saying that the limit in 
(1.3) exists. There still remain the problem of for which class of potentials V the 
expression (1.5) exists and whether (1.6) actually produces the propagator. 

We will demonstrate the usage of (1.5) and (1.6) on the <i-dimensional harmonic 
oscillator path integral. The harmonic oscillator plays a major role in quantum field 
theory and the recent advances due to Duru and Kleinert in the coulomb potential 
path integral(see [7] and references within). The harmonic oscillator carries the 

potential V (x) = m ^ x 2 , its internal functional integral is: 
(1.7) 



(- 



m 



exp 



. OJ+l 



m 
~2 



x 7 _i 



m 



■•" \2 2 



dx\ . . . dx u 



Notice that each d-dimensional integral factors into d products of one dimensional 
integrals, thous we shall compute the 1-dimensional harmonic oscillator internal 
functional integral. 

II. The Harmonic Oscillator. There are many ways to compute the harmonic 
oscillator path integral(see [7], [8], [9] and references within), we will use a method 
of computation similar to that of [8] . We differ from the popular techniques in that 
we will do the computation without prior knowledge of the classical path of the 
harmonic oscillator and we will rigorously do the computation with Nonstandard 
Analysis. As pointed out earlier, it turns out that properties of Classical Physics 
falls out naturally from a purely Quantum Mechanical derivation. 

7T 

It is well known that for < t < — , the propagator for the 1-dimensional 

A 

harmonic oscillator is 



(2.1) *(«,*.«) =(5^)' 



A 



sin At 



exp 



im 



A 



h sin At 



[(<7o 2 + q 2 ) cos At - 2qq ] 



Equation (2.1) carries a singularity at t 



TV 

A 



Singularity like that of (2.1) can 



in fact be given rigorous mathematical meaning if K(q,qo,t) is interpreted as a 
distribution(see [10]). Due to the form of (2.1), we would not expect the standard 



TV 



functional integral to exist at t = — as a function. 

A 

The popular method to compute the 1-dimensional harmonic oscillator Feyn- 
man path integral is by writing S [x (s)] = S [x cl (s) + Sx (s)] , where x cl (s) is the 
classical path of the harmonic oscillator which satisfies the equation of motions 
x cl (s) = —X 2 x cl (s), with the boundary condition x cl (0) = qo,x cl (t) = q. Namely, 



x cl (s) 



(2.2) 



q sin Xs + qo sin A (t — s) 
sin At 



The action integral becomes 



* m 

t 



(x cl ) 2 -x 2 (x cl y 



ds + 



m 
2~ 



{Sx) 2 - A 2 (SxY 



ds+ 



/ m (x cl Sx - X 2 x cl 5x) ds 
Jo 



Using Sx (0) = = Sx (t) and the equation of motions of x c \ the last inte- 
gral is after an integration by parts. Integrating by parts on the first integral 
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and using the equation of motions of x cZ , the first integral becomes — x cl x cl 

7Tl\ 

— \(q 2 + q%) cost — 2qq ] . Without much concern on the existence and mean- 

2 sin At L ' J 

ing of the path integral, we can write 



(2.3) 

/ ex P \ — k — r dx ( s ) = 

Jx(0) = q o I n J 

r i m \ *\ r5x(t)=o (iSlSxis)]} 



The path integral on the right-hand side is the quantum fluctuation; the integral is 
over all paths which starts from at s = 0, and end at at s = t. We leave it to 
the reader to look up the computation of the quantum fluctuation in the literature. 



We will give a rigorous treatment of (2.3) by using the time-sliced internal path 
integral in (1.5) and (1.6). In our work, we do not start with having knowledge of 
the classical path. We will start with an arbitrary bounded path w(s) which sat- 
isfies w (0) = qo,w (t) = q, and separate the propagator into a product of classical 
and quantum amplitudes. In this approach, we will see that the classical contribu- 
tion actually comes from Quantum Mechanics without prior knowledge of Classical 
Mechanics. 



To shorten the notation, we write 



(2.4) 

. n+l 



im \ 
We) 



m I Xj — Xj-i 
~2 



in 



n+l 



xl - 2x xx + x 2 n+1 - 2x n x n+1 + ^ 2x ? _ ^2 2x 3 x j-i ~ f2 A 2 ^2 



x~ 
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im \ , 
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1/ 
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Vo 













Vo 
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-1 








-1 
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-1 








-1 


2 -1 ... . 










1 
1/ 



+ 



-12-10 
0-120 
0/ 



>x = 



\ ( trp \ 



where T n is the (n + 2) by (n + 2) symmetric matrix, 

-1 ... 0\ 



(2.5) 



/ 1 

-1 


V 



S n 
-1 

-1 l-e 2 A 2 / 



with S n being the n by n symmetric matrix S n = A n — e 2 \ 2 B n , where 



(2.6) 



/ 2 -1 ... 
-12-10 
0-12-1 






V o 



\ 










0-12-1 
. .. -1 2 / 
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/I 







1 







B r 





Vo 



i 







1/ 



and x is the column vector 



(2.7) 



x = 



( x \ 

X! 



X n 



For notation convenience, we will use a bar instead of * to indicate the *-transforms 
of matrices, determinate of matrices, and vectors. 

We are interested in knowing when the internal functional integral of the har- 
monic oscillator exists. From (2.1), it would be reasonable to postulate the exis- 



7T 



tence of the functional integral for t < —. Indeed, this turns out to be the case. For 

A 

7T - 

t < — , ^(the *-transform of S n with n = u E *N — N) turns out to be *-positive 
A 

definite, which allows us to actually compute the integrals. 



Proposition 2.1. For < t < 



2 2 



\j A 2 (n + lf 



1 - 



7T 



12 (n + iy 



definite. 

Proof. An elementary computation shows that for k = 1,2, ... ,n, 



S n is positive 



(2.8) 



2/C7T 



sin 



n+l 



2k-K 



= 12 — 2 cos 



S7T 



n+l 



V sin St/ 



sin 



n+l 



Vsin^/ 



Hence, the n distinct eigenvalues of S n are 2 — 2 cos 

( ^fT ) - A 2 (£) - To show that 
S'n is positive definite, it is enough to find the values of t for which the eigenvalues 

are positive, or cos ^^pyj < 1 — - 2 . Since cos ^^jj < cos (^^rj f° r & = 

/ \ A 2 t 2 

l,2,...,n, it is enough to find t for which cos y^Tj < 1 — — j. By Taylor 
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expanding cos (^pr) about for the first 3 nonzero terms, we have 



(2.9) 
cos 



7T 



71+1 



2+2 



< 1 - 



t 2 < 



2 2 



A 2 (n + iy 



XH 
2ri* 



1 - 



1 - 



7T COS T] 



7T 



+ 



7T COS T] 



2(n + l) 2 4!(n + l) 4 



< 1 - 



\H 2 

2n 2 



12 (n + IT 



where < r] < 



7T 



n+1 



When t < 



n 2 7r 2 



\j X 2 (n + iy 



1 - 



7T 



12(n+iy 



, we have t 2 < 



7T 2 n 2 



A 2 (n+1)' 



7T COS ?7 

12 (n + if 



□ 



7T 



Theorem 2.2. Let t e t anrf < £ < — . For any u £ *N - N, S u is positive 
definite in the *-transformed sense. 

Proof. *-transforming Proposition 2.1 and setting n = u>, we have that is positive 
definite when 



(2.10) 



< t < 



co 2 n 2 



\ A 2 (u + iy 



1 - 



7T 



7T 

A 



+ h. 



n 



where h is infinitesimal. When t is standard and < t < — , (2.10) holds. □ 

A 



7T 



From here on, let us take < t < — . We will now proceed to separate the 

A 

functional integral into a classical part and a quantum fluctuation part. Suppose 
w (s) is an arbitrary path with \w(s) \ < oo for < s < t. Furthermore, let 

w (0) = qo, and w (t) = q. We make the substitution xj = w (^^+t j + Vj = w j + Vj 
(notice that yo = = y n +i since w (0) = xq = qo and w (t) = x n+ i = q). Using the 
fact that T n is symmetric, we have 

(2.11) 

x*T n x = (y + wf T n (y + w) = w l T n w + y t T n y + w l T n y + y l T n w = 

w l T n w + y t T n y + (T n w) 1 y + (w*T n y)* = w f T n w + y l T n y + (T n w)* y + ty*T n y = 

w*T n «; + y l T n y + 2 (T n iw)* y 

By using y = = y n+i and writing T n as 



(2.12) 



T — 



/ 


s. 
V o ... 













0/ 



+ 



/ 1 

-1 





V o 



-1 




-1 l-e 2 A 2 / 
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we obtain 

(2.13) x l T n x = w'TnW + y l T n y + 2 (T n wf y = w l T n w + y*S n y + 2p t y , 



where 



(2.14) 



Vn 

v o y 



, y 



P = S r , 



( Wl \ 

W 2 



Wn-1 



fvi\ 



\VnJ 



, W 



( W ° \ 



\W n+ i I 



( w \ 






\w n+1 ) 



S„w — w. 



and 



We then have the following: 



Lemma 2.3. Under the assumption on t, let e = — , then 

n 



(2.15) 



ihe/ 



2niht 



exp 



n+l 

E 

i=i L 



m [ xj Xj-i \ _ m 2 2 
2 j 



dx\ . . . dx 7 



exp 



^ (w*T n w - p'S^p) 



m \ 2 / 1 



2nihtJ V detS. 



Proof. Since S n is positive definite, it is invertible. Since S n is symmetric, the 
following holds 



(2.16) 



y*S n y + 2p t y = (y + S^p)* S n {y + S' 1 p) - p'S^p. 



Using our shorten notation in (2.13) and (2.16), the integrals in (2.15) is equivalent 
to 



(2.17) exp 



im 
he 



(w*T n W - p'S^p) 



m 



) 6X ^ V ~2h~ Z I 1 " z 



2irihe 



In obtaining (2.17), we performed the change of variables Xj = Wj + yj, and then 
from yj + (S~ 1 p)^. = Zj. We get (2.15) after diagonalizing S n and doing the 
decoupled integrals. □ 

Corollary 2.4. Under the previous definition ofw(s), The 1- dimensional har- 
monic oscillator internal functional integral is well defined and it is equal to 



(2.18) 



exp 



im 



/ m \2 / 1 
\2nihe) V detS,, 
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where the bars denote -^--transform. 

Proof. This is just the *-transform lemma 2.3. □ 

Remark 2.1. There is no restriction on the choice of the path w except that 

it starts at qo and ends at q. We will show that the exponential part of (2.18) 

turns out to be the classical amplitude of the previous formal calculation of the 

propagator and the other factor is the quantum fluctuation. If we choose w to 

be the *-transform of the classical path, it can be shown that each entry in p is 

iv^T vu A 

infinitesimal, and — is infinitesimally close to — \(qo 2 + q 2 ) cos Xt — 2qqo] . 

e sin At LV ' 

There are many techniques to compute the quantum fluctuation lim \ — — — 

n^oo y edetbn 

in the literature; we present a rigorous method to compute the limit with Nonstan- 



dard Analysis by computing st 



1 



edetS u 



for uj e *N - N. 



Proposition 2.5. Let e 

j matrices given by 



t 



n 



Denote A^ n and Cj >n with < j < n to be the j by 



( 



(2.19) 



A ■ — 



-1 



-1 






V o 







1 



\ 






•• 

2 -1 
-1 1 J 

. 0\ 
. 



, and 





Vo 



1 
0/ 



Define Dj n = det \\Aj :n — e A C n . After * -transforming, we have that for to E 

ktX 

*N - N, D KuJ w *cos( ) = *cos(fceA). 



Proof. For k = 1, Ai jU 
(2.20) A 2jU , = 



(1), C l uJ = (0) and D ljU = 1. For k = 2, 
2 -1 



-1 1 



G 



2,w 



1 





and D 2)LU = 1 — e 2 A 2 . Hence, the claim is true for k= 1,2. We expand Dj jn on the 
top row, and get the recursion relation 



(2.21) 



Dj,n = (2 - e 2 A 2 ) D 3 . x 



D 



j-2,n, 



2<j<n 
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Since we are interested in w e *N — N, we will consider the cases for which 4 — 
(2 — A 2 e 2 ) > 0. To solve the difference equation, we substitute Dj n = Aa^ 1 into 

(2.21) and get a 2 - (2 - A 2 e 2 )a + 1 = with solutions 

(2 - A 2 e 2 ) ± i\U- (2- A 2 e 2 ) 2 

(2.22) a ± = ^ 1 \ ^ '-. 



Both solutions a± have norm 1. Thus, we can denote a± = e i±e , where 9 = 
arg (a+) > 0. Hence, 

(2.23) Dj, n = A+exp {i (j -1)9} + A~exp {-i (j -1)9}. 
Solving for the initial conditions 

(2.24) D 1 = 1 = A+ + A~ 

D 2 = l- e 2 A 2 = A+a+ + A~a-, 

we get 

(2.25) A - : i ( XC 



.2V4- A 2 e 2 / ' 

We now proceed to get an estimate for 9. By definition of 9, cos 9 = 1 — 



2^2 



A 2 e 



2 

02 ^2 £ 2 

After expanding the left-hand side about 0, we get 1 — cos r] =1 — , < 

rj < 9. To estimate 9 by Ae, we write 9 = Ae + 0, and obtain 

(2.26) i-J^cc^l-^, => 

= — Ae ± ~ - =^ </> = — Ae + 



^/cos 77 ' ^/cos ?y ' 

where the last implication is due to 9 > 0. Thus, 9 = Ae — Ae I 1 



^/cos 77 / 

By *-transforming the above and setting j = k, n = u, 2 < k < u>, u e *N — N, 
we get 

(2.27) e = - « 0, 9, r\ « 0, 1 = « 0, « A" « -. 

to V *cos rj 2 

Equation (2.27) and (2.23) implies that 

(2.28) D kjLJ « *cos{(£;- 1)0} = 

*cos j (^Ae - - fcAe ^1 - ^== ) ) } ~ * cos (^ Ae ) • 

In the last ~ in (2.28), we used the fact that the cosine function is uniformly 
continuous, which translates into *cosx ~ *cosy whenever x ~ y in the language 
of Nonstandard Analysis. □ 

With the aid of Proposition 2.5, we can show the following: 
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sin (At) 



Theorem 2.6. edetS^ is infinitesimally close to 



A 



Proof. As before, we will use bars to denote the *-transform of determinate of 
matrices. Let Sj yU = det ||A,-, n — e 2 A 2 5j iTJ ,||, where Aj yU is j by j matrix as defined 

in (2.19), Bj n is the j by j identity matrix, e = — , and 1 < j < n. Notice 



n 



that S njU = det\\S n \\. Expanding Sj iU and Dj n on the bottom row, we get the 
recursion relation Sj jn = Dj jU + (l — e 2 A 2 ) Sj-i tn , or equivalently, Sj tH — Sj-\, n = 
Dj, n — e 2 A 2 Sj_i )n . Summing the last equality gives 



(2.29) 



S n ,n Sl,n — Dj,n e ^ ^ ] Sj — l,ni 
3=2 j=2 



3=2 3=2 

From the recursion relation, we also get that for j > 3, |5j_i jn | < ^X)fc=2 \Dk,n\) + 

\Sl,n\- 

We now *-transform the second equation in (2.29) and write from proposition 
(2.5) t) m ,uj = *cos (meA) + h m , where h m is infinitesimal. We get 
(2.30) 

edet S u = eSi^ + ^ e*cos ( V^h. j _ e * C os (eA) + e ^ h m — e 2 A 2 ^ e ^m-i,u;- 

m=l \ u / m=2 m= 2 

The set {/i m |l < m < uj} is an internal set, so max h m G {/i m |l < m < u>} and it 

l<m<w 

is an infinitesimal. Thus, e^ =2 h m < ( — J oj max /i m 0. From the bound on 

\0J J l<m<o> 

\Sj- 1;n \, we get 



(2.31) 



|e 2 A 2 £ c5 m _ 1>w | < e 3 A 2 |^| + e 2 A 2 £ e 



m=2 

e 3 A 2 |^|+e 2 A 2 















+ \Si, u \ 


h 


Lm=3 









I> +f|Si, 



\m=2 



0. 



Finally, using the limit of Riemann sums in the language of Nonstandard Analysis, 

we have Ylm=i e * cos i 11 ^) ~ Jo cos ( / ^' s ) ^ s = Since the other 2 terms in 

A 

(2.30) are also infinitesimals, the result follows. □ 

We are now ready to derive the classical amplitude from the exponential in (2.18) 
by using results from theorem 2.6. 

Proposition 2.7. The exponential in (2.18) satisfies the following 

im A 



(2.32) exp 



exp 



h sin At 



[{ql + q 2 ) cos \t - 2qq ] \ , 
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where qo = xo, and q = x^+i 

Proof. From the definition of p,w,T n and S 1 " 1 in (2.12) — (2.14), we multiply out 
p t S~ 1 p and express it in terms of w t T n w. 

(2.32) 

p'S-'p = (w'Sn - w ) S- 1 (S n w -wj = (w* - wS- 1 ^) (S n w - vj 

v^SnW — "ufw — w w + w S^ 1 ')!; = w t T n w — Wq — (l — e 2 A 2 ) w^ +1 + 
(Sn 1 )^ + w w n+1 (S- 1 ) ln + w w n+1 (5" 1 ) nl + w 2 +1 (^ 1 ) nn 

After *-transforming (2.32), we get 

(2.33) 

1 (w*f u iv - p'S-'p) = \wl + (1 - e 2 A 2 ) w^ +1 - wl (S- 1 )^ - 



1 



+i {S u 1 ) uu -w w u+1 (S u 1 ) lu -iv iB u+1 {Sj) ul ] = 

<ll ( 1 " + ? 2 ( 1 " " 2 gg0 ( " ir+ , 1( " irl " eW 

2 2 _ 



edetS u 
A 

sin At 



- [<Zo — ^-l,^) + <? 2 (<S^,^ — S^-i^) — 2gg ] — eA 2 g 

[(<?o + I 2 ) cosXt - 2 <?<?o] • 



In the third line above, we used *-Cramer's rule and our previous definition of Sk,u- 

The factor (— 1) <JJ+1 comes from a cofactor expansion in the numerator of Cramer's 

rule for (S^ 1 )^ and (<S'j 1 ) a , 1 ; the factor (— l)^ -1 comes from the determinate of 

a triangular matrix with — l's along the diagonal after the latter expansion. In 

the fourth line, we used results from Theorem 2; namely, <S WjW = detS^, edetS^ ~ 

sin At - - - > r-i 

— - — , and S u ,u - <Sw-i,w ~ A^,w ~ cos At. □ 
A 

Notice that in (2.33), the end result dependent only on the endpoints of the path 
w. It does not matter which path is choosen as long as it starts at q and ends at 
qo. Hence, it is not necessary to use the classical path x cl to do the computation. 

Theorem 2.8. Fort < —, the one dimensional harmonic oscillator standard func- 

X 

tional integral is given by 



Proof. This follows from corollary 2.4, theorem 2.6, and proposition 2.7. □ 
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Corollary 2.9. For the d- dimensional harmonic oscillator standard functional in- 
tegral, we have 

(2.35) 

St \ J \2nihe) 




Proof. Follows from factoring (2.35) into products of one dimensional harmonic 
oscillator standard functional integrals and theorem 2.8. □ 
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